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We investigate quantum tunneling in a translation invariant chain of particles. The particles 
interact harmonically with their nearest neighbors, except for one bond, which is anharmonic. It 
is described by a symmetric double well potential. In the first step, we show how the anharmonic 
coordinate can be separated from the normal modes. This yields a Lagrangian which has been 
used to study quantum dissipation. Elimination of the normal modes leads to a nonlocal action of 
Caldeira-Leggett type. If the anharmonic bond defect is in the bulk, one arrives at Ohmic damping, 
i.e. there is a transition of a delocalized bond state to a localized one if the elastic constant exceeds 
' a critical value Ccrit- The latter depends on the masses of the bond defect. Superohmic damping 

occurs if the bond defect is in the site M at a finite distance from one of the chain ends. If the 
observation time T is smaller than a characteristic time tm ~ M, depending on the location M of 
OO ' the defect, the behavior is similar to the bulk situation. However, for T 3> tm tunneling is never 

CN| ' suppressed. 
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PACS numbers: 03.65.Xp, 05.30.-d, 61.72.-y 



I. INTRODUCTION 



' The influence of environmental degrees of freedom (DOF) on quantum phenomena, like e.g. tunneling, has been 
of great interest during the last decades [H, 0, 0] ■ ^ significant progress came from the investigation of a more or 
I less phcnomenological model where a particle in a one-dimensional potential V{q) or a free particle is coupled to a 
bath of harmonic oscillators with the spectral density JU^o)- The quantum dissipation generated by the bath depends 
qualitatively on the low frequency behavior of ./(w) [H, H, Q. A particular interesting case is Ohmic damping, when 
^ I J{uj) ~ w for low enough frequencies. In that case and for a symmetric double well potential, the particle at zero 
temperature undergoes a transition from a delocalized state to a localized one if the coupling constant between the 
particle and bath exceeds a critical value [3| . An interesting observation has been made by Caldeira and Leggett 0] . 
The exponent of the exponential factor for the tunneling probability is multiplied by 77, the phcnomenological friction 
' coefficient of the corresponding classical dynamics. This relationship between classical and quantum dissipation 
has been deepened and generalized by Leggett for an arbitrary linear coupling between the particle and bath 
10 ■ coordinates. 

\ Let q{uj) be the Fourier transform of the classical particle trajectory q{t) and 
O ' ~ dV 

i^oHg» + 7pH = o (1) 

O ■ 

the transformed classical equation of motion, where Kq{uj) contains the dissipative influence of the bath. Then 
k> the reduced Euclidean particle propagator G£;(q', T|q, 0) (where the harmonic DOF have been eliminated) can be 
, represented by a path integral in the imaginary time t = —it [3] 

^. 

GE{q',T\q,0)^ I I?[q(T)]e-*^['Ml . (2) 



q(0) = q 
<l{T) = q' 



The action 

S[q{T)] = So[q{T)] + Snonlocaliqir)] (3a) 



contains the local 



So[q{T)] = / dr 



^q{rf + V{q{r)) + gg(r)^ 



(3b) 
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and the nonlocal part 

T r 

Snonlocal[q{T)] = - J dr J dr' K {t - T')q{T)q{T' ) = -i-i / dio K(u;)\q{u;)\^ . (3c) 


The second equality holds for T — oo. Kiuj) is the Fourier transform of the integral kernel K{t) and it is related to 
K[){u}) by K{lu) = ^KQ{—i\ui\). If the Euclidean Lagrangian of the particle-bath system is 

L = Li) + Li , Li = Lbath + Lint (4a) 

with 

Lo{q,q) = \Mpe + V{q), (4b) 



and 

jv r , ^ ^ 2' 



Li{q,x;q,i) = | X! 

then0,i,[3] 



N 

Too 

a=l 



(4c) 



-Id^^ (5a) 







and 



TT 



,co sh(^(f -M )) 
sinh(a;Y) 



K{t) = - I dtoJiiu) . y (5b) 



with the spectral density 



N 2 



a=l 



For a finite T the kernel K{t) and the paths ^(t) can be periodically continued. Then the Fourier scries for K{t) is 
given by the Fourier coefficients Q 

Kn-^Y.—^^ (6) 

with Vn = (27r/T)n, n = 0, ±1, ±2, .... Here Mp and {toq} are the masses of the particle and harmonic oscillators, 
respectively. {uja\ are the oscillator frequencies and {cq} are the coupling constants between q and the coordinates 
of the oscillators {cCq,}. Note that {xa} are not necessarily positions, but can represent normal mode coordinates of 
vibrations, etc. 

In the following we will restrict ourselves to a system of N particles whose potential energy V{xi, . . . , afjv) includes 
harmonic and anharmonic interactions. Without an external field, V must be translationally invariant. However, 
since the coordinates of Langrangian (4) are not specified, it is not necessarily invariant under translations. This 
has motivated Chudnovsky [g| to apply the Caldeira-Lcggett approach to a system of two particles (i = 1,2) with 
positions Xi and masses Mi, coupled to oscillators with coordinates Xa, frequencies Wq, and masses TOq,, a = 1, . . . , N . 
The corresponding Euclidean Lagrangian is of the form 



AI AI 1 ^ 

L = -^xj + -^x\ + V{xx - Xa) + - E TOa {x\ + w^(xc, - xif\ (7) 



with V {x\ — X2) being the interaction energy between both particles. The coupling of particle i = 1 to the oscillators is 
assumed to be zero. It is obvious that L is translationally invariant under provided that 
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{xa} are real space coordinates. Surprisingly the elimination of harmonic DOF does not lead for zero temperature to 
a nonlocal action of Caldeira-Leggett type (cf. Eqs. ([3c]), (j5bp .([5|)). The Fourier coefficients of Chudnovsky's kernel 



have the form 



k;, . 1 MM ^, (8, 

Ml +M2+y ° " 



which differs quahtatively from the Caldeira-Leggett type result ([6]). However, for Mi oo one gets 



^ I (A/i ~ Af2)z.^ -l-Y. 2"i«^a (9) 



where the last term, up to the additional factor — t'^, corresponds to the second order time derivative of the kernel, 
identical to A'„ from Eq. ^ , if one chooses = \plmaiOa [1| • 

Then the question arises: Does a translationally invariant model in general lead to a nonlocal action, which is not 
of Caldeira-Leggett type? This is one of the main points we want to investigate among others in our paper. It will be 
done for a microscopic(within a Born-Oppenheimer approximation), explicitly translationally invariant lattice model 
with a defect which cannot diffuse. We will show how the normal mode coordinates for the harmonic DOF can exactly 
be separated from the anharmonic ones. This leads to a Lagrangian of the form of Eq. (4) with coupling constants 
Cq, frequencies LOa and a spectral density determined by the microscopic model parameters. Such a microscopic 
justification of Lagrangian (4) has been presented for quantum diffusion @. There a particle diffusing through an 
elastic lattice is considered. If, however, that particle cannot diffuse and is an integral part of the lattice, e.g. an 
impurity which can tunnel only between two positions, one has to separate the center of mass (COM) and relative 
coordinates of all particles. For such a situation, Sethna [To| has estimated the coupling constants Cq, by comparing 
the strain field of an elastic monopole of an impurity with the displacement for a longitudinal mode. But, as far as 
we know, there is no microscopic derivation available for the quantities Cq,, Wq, and J{lo) for a non-diffusing impurity 
in a lattice. Besides such a microscopic derivation we will show that the quantum behavior of the defect is sensitive 
to its location. 

The outline of our paper is as follows: Section [TT] presents our model and outlines the main steps leading to 
Lagrangian (4). The implications for the quantum behavior will be discussed in the third section with a special 
emphasis on the role of defect location. A summary and conclusions are contained in the final section HVl Appendices 
A, B and C contain details on the separation of the harmonic and anharmonic DOF. 



II. MODEL 



We consider an open chain of N particles with masses m„, n = 1, . . . , iV and harmonic nearest neighbor interactions 
and one anharmonic bond representing a defect. This model could describe a linear macromolecule with an impurity. 
The classical Hamiltonian reads 

^ 1 

H = J2o pI + V{xi,...,xn) (10a) 

^ 2m„ 

n— 1 

with the potential energy 

C ^"^ 

V{xi, ...,xn) ^ — y {Xn+l ~ Xn - a)'^ + Vo{xm+i ~ xm) ■ (10b) 

71=1 

in which Xn is the position of n-tli particle, C is the elastic constant of the harmonic nearest neighbor interaction, 
a is the equilibrium length of the harmonic bonds and Vo(a;A/+i ~ xm) is the energy of the anharmonic bond defect 
located between the sites M and M + 1. The potential energy V is explicitly translationnally invariant. Several ways 
exist to separate the harmonic and anharmonic DOF, which finally lead to the Langrangian (4) (see Appendices [Bl 
and lCj) . One (see Appendix IB|) . which is also applicable in higher dimensions, is to introduce the COM 



Xd = ' {rriMXM + niM+iXM+i) 



(11a) 
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and relative coordinates of the bond defect 

QM = XM+i - XM ■ (lib) 

Then, xi, . . . , xm-i, X^, x^z+i, . . . ,xn arc harmonic coordinates. Their kinetic and potential energy can be diago- 
nalized by introducing normal coordinates, qm is linearly coupled to these coordinates. As a result, one obtains the 
Lagrangian (4). Here we will choose a different approach applicable to Id systems, which, however, leads to the same 
Lagrangian. Let 

^ N N 

Xc = — ^ ninXn , Afc = ^ m„ (12a) 

^ n—l n—1 

be the COM of all the particles and 

qi = Xi+i - Xi - ai , z = 1,...,7V- 1, 

o-M = , and ai ~ a otherwise, (12b) 
be the relative coordinates, respectively. Using notation qo ~ X^ Eq. (|12ap . (jl2bp has the form 

N 

qi + at ^ AinXn , i = 0, 1, - 1 . (13a) 

n=l 

with 

Atn = ^'^O.i + iSi,n + l - <5j,„)(l - So,i) 

Let TTi be the canonical conjugate momenta of qi. It is easy to prove that Eq. (|13ap implies 

JV-l 

Pn = ^ AinTTi , ?1 = 1, ...,N . (13b) 

1=0 

Substituting pn into Eq. (|13a[) yields 

1 2 1 1 1 \ , 1 

^ 2M:"° + 2 ^ + ^ j " it 
C 

+ 2 ^ + • (14) 

(i#M) 

Here the first term is the kinetic energy of COM, which will be dropped from now on. Note that the use of relative 
coordinates introduces a coupling between the momenta. In the next step, we perform the canonical transformation 

T^i^Pi, i = 1, ...,iV - 1; j ^ A/, 

TTM = Pm H — {niM+ipM-i + niMpM+i) , (15a) 

rriM + niM+i 



q, = qt, i^l,...,N-l;i^M±l 



qM-i = qM-i ■ qM, (15b) 

Tom + tom+1 



rriM 

qM+i — qM+i ; qM 

VIM + TOA/+1 
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This leads to 

H = Hd + Hharm + Hint (16a) 

where 

2mMmM+i 2{mM+mM+ir 

is the defect Hamiltonian, 

N-l „ N-l 



Hd=-z Pa/ + ^o(gM) + TTT ^ (16b) 



1 C 

Hharm = 9 X! '^iJ^iPj + ^ XI " 

is the harmonic part of the Hamiltonian and 



2 ^ ^'J-'-J 2 

i,j=i 1=1 



i?mt = -C ■ qM (16d) 

is the coupling between the two. 

The matrix T = {Tij) in (|16cp depends on the masses rrii (see Appendix Let e^"^ = 

\ . . . , e^^^-j^, e^"^j^, . . . , e^l^)* and Aq, be the normalized eigenvectors and eigenvalues of T, respectively, and 

S = (e(i),...,e(^-2)) 
be the orthogonal matrix, which diagonalizes T, i.e. 

S-iTS = A, Ko.p = K5c.p- (17) 
Then we can introduce the normal coordinates 

N-l N-l 
Xa = ^ Saiqi , Pa = ^ SaiPi (18) 

i=l i=l 
i^M) (MM) 

for a = 1,2, N — 2 such that Hharm becomes diagonal 

2 



N-2 

Hharm ^^Y. t^"^' + ' (1^^) 



a=l 

The interaction term takes the form 

N-2 



Hint = - X CaXaqM ■ (19b) 



a = l 

The M dependent coupling constants are given by 



Ca — C- 



(") , (a) 



(20) 



rUM + TOA/+1 

The final step is the Legendre transformation of H , which leads to the Euclidean Lagrangian 

L = Ld + Li , Li = Lharm + Lint (21a) 

where 

Ld = ^mqlf + ^o('?m), (21b) 

and 

N-2 

2 



^ N-2 

Li = i^J^^'^ 



a=l 



^ 2 

2 I , ,2 ' 



+ UJ^ \ Xa ~ 2 qM 



Co 



(21c) 
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Here we have used the equahty g^/ = Qm (cf- Eq. (jl5bp ) and 



N-2 



r,2 



y-^^C^I^^i^^^, (22) 



which follows from the completeness of the set e*^"' of eigenvectors and l/(mQ,(jj^) = 1/C (see below). Eq 
allows us to include the countcrtcrm — —. — q\, in Eq. (I16bl) into Li. This counterterm, the role of which 
has been discussed by Caldeira and Leggett [ll|, results from the canonical transformation, Eq. (|15ap . (|15bp . This 
transformation eliminates the coupling between the momenta of the harmonic DOF and that of the defect and generates 
coupling between the normal mode coordinates {xa} and the corresponding defect variable (cf. Eq. (jl6dp ). Due to 
the disappearance of the counterterm, there is no frequency renormalization for the bond defect 0]. The Lagrangian 
(21) is identical to that of Eq. (4). The masses TOq and frequencies u)a follow from 

ma = -r— , = [CXa)^ (23a) 



and the reduced defect mass to is given by 

toa/tom+i 



TO A/ + TOM+l 



(23b) 



Note that these results are exact for one dimensional systems. It can be shown that for an arbitrary impurity in a 
two- or three-dimensional system the Langrangian (21) can be derived within a kind of harmonic approximation [l^ . 



III. QUANTUM TUNNELING 



In this section we will investigate the zero temperature quantum behavior of the anharmonic bond defect embedded 
in the harmonic chain as shown in Fig. [TJ We will assume that Vo((7j\/) is a symmetric double- well potential with 



• • • # # • • • • (I) 

1 M M+1 N 

a ai 

• • • # % • • • • (11) 

1 M M+1 N 



FIG. 1: Two degenerate classical ground states of the open chain with A'^ particles. The masses m„, n 7^ M, M are chosen to 
be equal, a is the equilibrium length of the harmonic bonds and 03, ai the two degenerate equilibrium lengths of the anharmonic 
bond. 



degenerate minima at q^j = as > and q^^ = ag > ag. Then the classical ground state of V{xi, . . . ,xn) (cf. Eq. I10b[) 
is twofold degenerate (see Fig.l). Therefore the low-lying eigenstatcs form doublets. Neglecting the excited doublets 
at zero temperature is justified if the bare tunneling splitting of the ground state doublet is much less than the 
frequency of the upper phonon band edge luq (see Eq. (HSUDd. If the total number N of particles in the chain 
is macroscopically large and the particle number M is in the bulk of the chain, i.e. M = 0{N), one might have 
expected a suppression of tunneling since the change from, e.g. q'^j = to q'^j = ai, would require a translation of 
the macroscopic mass of the left and right harmonic parts of the chain. We will see that this naive expectation is not 
always correct. 

On the other hand, if the defect is close to one of the free boundaries, i.e. either M = 0(1) or (iV — M) = 0(1), 
only a finite mass ~ M has to be translated. Consequently tunneling cannot be suppressed. This qualitative M- 
dependence should follow from that of the kernel K{t), which itself results from the strong Af -sensitivity of the 
spectral density J{uj). Subsection IIII A] discusses this phenomenon. The limit toa/ 00 (or toa/+i — > 00), motivated 
by the conclusions drawn in Ref. IS], will be discussed in Subsection IIII Bl 
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A. Location-dependent quantum tunneling 

As we will see in Subsection llll BI tlie form of the nonlocal action Eq. pc]) does not depend qualitatively on the masses 
TO„, even if one of the defect masses niM or iriM+i becomes infinitely large. Therefore we will choose for convenience 
TO„ = m. Let us start with the situation where the bond defect is located within the bulk, i.e. M = 0{N), so that 

O < lim ^ =(,<!. (24) 

One can prove that the tunneling phenomena do not depend on ^ if it is different from and 1. Therefore we choose 
M = N/2 and without loss of generality to be even. This choice and the assumption m„ = m allows us to determine 
the eigcnfrequencies and the eigenvectors e'"\ a ~ 1, . . . , iV— 2 exactly for finite N. A calculation, whose technical 
details are presented in Appendix [XI results in 



+(«) - 

^ " V A - 1 



( ^ ) ^ = '^osin ( yj , wo = 2y— , (25a) 



'sm{q+j) , 1<J<M-I = f-1, 

sin(<7+(j-l)) , M + l = f + l<j<Ar-l, 



e/"^ = a/| sin(g-j) , l<j<N-l;j^M^^, (25b) 



where 



TT 



-(2a -1) 



A- 1 

= { (25c) 

TT 

for a = 1, . . . , A/2 — 1. It is easy to see that e^^"'' are the symmetric eigenvalues with respect to j —t N — j and e"*-"-* 
are the antisymmetric ones. From this, m„ = m and Eq. (|20p it is obvious that the bond defect does not couple to 
the antisymmetric vibrational modes, as may be expected from the symmetry of the problem. With these results we 
can calculate the spectral density. Making use of Eqs. ((20)) . (|23a[) . (|25ap . (j25bp and (|25cp we get from Eq. ([5c|) in the 
thermodynamic limit N oo that 

TT 

J (^) = ^C'^o j dq cos^ (^1^ sin (^^^ d(uj -ivq sin 




(26a) 



In the limit w <C we obviously have 

J{uj) ^ C— , (26b) 

which corresponds to the situation of Ohmic damping. The Ohmic damping results from two facts: First, the density 
of states 5(w) of the vibrational modes in a one-dimensional lattice is constant for uj <^ loq and second the squared 
couphng constant involves the factor sin^(gQ,M) which, for A^ ^ oo, Af ^ cxd with M/N = 0, 1), oscillates faster 
and faster so that it can be replaced by 1/2. It is emphasized that the condition M = 0(1), i.e. ^ = or 1, will 
change the shape of J{lv) qualitatively. 

Now we can calculate the kernel K{t). For looT ^ 1 and |wo(5 ~ kDI ^ 1 the fraction in Eq. (|5bp can be well 
approximated by exp(— wo|t| sin |). The influence of the oscillators on tunneling of qm is determined by the large 
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— T behavior, i.e. by the low frequency modes. Therefore, substituting J{oj) from Eq. (j26bp into Eq. (j5bp we find, of 
course, the well-known result for Ohmic damping [l], [3, 

K{t) ^ -CcJo , ^ ,o , WOT > 1. (27) 

As a consequence, there exists a critical elastic constant Ccrit so that the anharmonic bond can tunnel for C < Ccrit, 
despite macroscopic masses have to be moved (see Figure 1). For C > Ccrit symmetry is broken. If the anharmonic 
bond is prepared in its ground state, e.g. = a^, it will remain there on average. 

If, however, the bond defect is located close to one of the boundaries, so that either M = 0(1) or N — M ^ 0(1) 
(^ = or 1 in ([24])), the situation changes. Note that we perform first the thermodynamic limit ^ oo. Then 
M = 0(1) means that M may equal 1, 2, 10^ or even a larger, but still finite number. For to„ = m one can easily 
show that 

1 „ / (a) , (a) 



-C [eiy_, + e\y^,) = CK. sin (g„M) (28) 

where Na is the normalization constant of {ei"^} (see Appendix [A|) . Replacing Ma by its low frequency behavior 
(2/A^)i/2 and taking the limits loqT > 1 and \ujq(^ - \t\)\ > 1 yields 

Km{t)^]-Cuj^- f dqqsm^iqM)e-^'^"^^^^ (29) 
2 n J 


for the kernel. In order to be consistent we also replaced uj{q) = ujosm{q/2) by its low frequency dispersion uj{q) = 
^ujoq- The integrand of Eq. ([29]) involves two g-scales 

1 1 

M ujo\t\ 

Equating qM = Qt defines the timescale 

TM - cu-^M. (31) 

The physical meaning of tm is as follows: The path integral formalism 0, [3 allows one to investigate quantum 
tunneling by determining the instanton solutions, i.e. the solutions of the classical equation of motion for a double 
well potential in imaginary time. The width r^-mk of a single instanton is 



^kink 



1/2 



If we assume that ^^"(9^/) ~ C, the elastic constant of the harmonic bonds, then Tkink ~ so that 

Tm ~ MTkink ■ (32) 
The r-dependence of Km{t) is sensitive to whether \t\/ Tkink ^ M or vice versa. Let us start with the long time limit 
(i) \t\/tm > 1. 



Then it follows from Eqs. ([30]) and ([3T]) that qr "C qM- The major contribution to the integral in Eq. (|29|) comes 
from q < q-r qu- Therefore, we are allowed to replace siT?[qM) ~ s\T?[q/qM) by {q/qM)^ = ilM)"^, which 
leads to the spectral density J{uj) ~ uj^ at low frequencies corresponding to superohmic damping. It implies 
that Km{t) ~ T~*. The precise result is 

48 1 

KM{T) = —Cu:a- -M\ |T|»rM . (33) 
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1 10 100 1000 10000 



FIG. 2: r-dependence of Km{t) for A4 = 5, 10, 40 and 160 (from bottom to top) on a log- log representation. The dotted and 
dashed line, corresponds to r^^ and r^'', respectively. The crossover at r « tm from the behavior to that of t~* can 
clearly be seen. The parameters have been chosen as foUows: T — 10^, C — m— I^luq — 2. 



(ii) A different situation takes place for |t|/tj\/ <C 1. 

Then q^- 3> Qm and the integral in Eq. (|29p must be decomposed in two contributions J dq ■ ■ ■ + J dq ■ ■ ■ . The 

qM 

luqt/M. For the second one we are allowed to 
It ^ 9m so that the function 
exp (— i g/g,- ) varies 

slowly. Accordingly t ^ tm corresponds to an 'effective' spectral density J{uj) « uj, i.e. to Ohmic damping. It 
is straightforward to estimate both integrals. As the final result we obtain 



first integral yields a constant in the leading order in ojqt/tm 
replace siT?{qM) = siv?{q/qM) by 1/2 since the main contribution comes from q 
siT?{qM) = siv?{q/qM) is oscillating fast between zero and one, whereas exp (— ^wolrlg 



Km{t) = -Cujq 

TT 



1 



8M2 



l + O 



( 



\ M 



1 



(c^or)2 



(34) 



for Tfcinfc ^ T ^ TM ~ MTkink- Taking the limit AI ^ oo in Eq. ((34|) restores the result ([27|) for the kernel 
K{r). 

For M = 0{1) or N — M = 0(1), and TV ^ oo we obtain for Km{t) a crossover at tm from the power law 
for T <C Tm to r^^ for t ^ tm- Figure 2 illustrates this behavior for Km{t), calculated numerically. On the log- log 
plot of Figure 2 the crossover between both power laws can easily be observed. 

This crossover is related to the M-depcndence of the spectral density because the defect-phonon coupling constants 
Ca are M-dependent. If the 'observation time' T (which is the time T in Ge from Eq. ^) is smaller than tm then 
(r — r') will be smaller than tm, as well. Consequently the kernel Km{t — t') entering the nonlocal action (Eq. ((3c|) ) 
decays as l/(r — r')^. If, however, T is larger than tm then it is possible that (r — r') becomes larger than tm, as 
well. For those values of (r — t') the kernel Km decays as 1/(t — r')^. This discussion reveals that the choice of the 
'observation time' T allows to fix the 'large'-r behavior of Km{t). where, of course, t <T. If M is far away from the 
chain end the crossover time tm is correspondingly large. Increasing M even more tm increases, too. Nevertheless, 
the 'observation time' dependence still exists. It disappears for M ^ oo, only. We remind the reader that the limit 
N ^ CO has to be taken first. 

Now making use of the analogy P, [^, 0] between the calculation of the action of a multi-instanton configuration 
interacting via Km{t — t') and a one-dimensional Ising model with the coupling constants Jnm ^ KM^n — m) ^ 
\n — m|~*' we can conclude the following: If the 'observation time' T is smaller than tm then we have Km{t) ^ t~^, 
i.e. Ohmic damping. In that case the bond defect may tunnel for C < Ccrit{T), whereas symmetry becomes broken 
for C > Ccrit{T). Note, that this is not a sharp transition at Ccrit{T) since finite T corresponds to a finite Ising chain 
which does not exhibit a sharp phase transition. What really happens when increasing the coupling constant C is an 
increase of the correlation 'length' ^(C). As soon as ^(C) equals the 'size' T of the Ising-chain a 'long range' order 
occurs. However, if T is much larger than tj\/ « MTkmk we have Ki^i(t) ^ t^^ and tunneling is never suppressed 

One might be puzzled by these conclusions since the transition for Ohmic damping to decoherence for M = 
0{N),N oo occurs for large T, or to be more precise it becomes a sharp transition for T = oo, only. As we 
already stressed above, the transition for T < tm is not sharp. The relevant phonons contributing to Km{t) have 
wavenumbers 5 ~ <?,- ~ <?t = 7^- This makes the 'effective' spectral density Ohmic. Mapping the situation for 
M = 0(1) again onto the Ising chain of length L results in the coupling constants Jnm decaying like Jo\n — m|~^ 
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for |n — m| < M and like Jo|n — m\~^ for |n ^ to| > M . It is clear that there is no sharp phase transition for finite 
L. But if L < Af , the coupHng constants decay as Jo|n — m|~^. For a fixed temperature T (not to be confused with 
the 'observation time' T) there will be no magnetic order (= coherent tunneling) for Jq < JQ,crit{T){^C < Ccrit{T)). 
For Jo > Jo,crit{T){= C > Ccrit{T)) a crossover to the 'long range order' (= decoherent tunneling) takes place. 
Accordingly the quantum tunneling phenomenon is richer for M = 0(1) than for M = 0{N). 

Actually we may think also in the real time t terms that as long as only the phonons with relatively high frequencies 
w ~ 1/i and hence short wave length such that 1/q < M participate in the interaction with the anharmonic defect, 
the latter does not 'feel' that the chain is finite and behaves as in the Ohmic case. In the course of time the lower 
frequency phonons with higher wave length 'reach' the end of the chain and a crossover to a superohmic behavior 
takes place. 



B. Dependence on the masses of defect 



In this subsection we will assume that 



n ^ M,M + 1, 

n = M, 

n = M +1. 



(35) 



For A-I = 0{N) we may choose without loss of generality M = N/2 and N being even. It is easy to prove that e^''"'' 
from Eq. (j25b[) remain eigenvectors of T with q^, a — 1,2,..., N/2 — 1 given by Eq. (|25cp . The remaining {N/2 — 1) 
eigenvectors are of the form 



sin(gaj) 



, l<J<f-l, 



(36a) 



smiq-{N~j)) , f + l<j<iV-l 



Qq, is a solution of transcendental equation. Let us introduce the quantities = m/Mi, i ~ 1,2 and S = /3i/?2/(/3i+/32)- 
The limit Mi —^ oo (or M2 00) implies S 0. Since the results in [8] motivate us to study, e.g. M2 —* 00, we find 



1 



26 



la = la 



N -1 



tan(^) 

in the limit 6^0. Substituting Eqs. ([Ma]) . ((36bl) into Eq. ^ yields 



0(6^) 



(36b) 



in the leading order in S. As a result J{uj) and therefore K{t), too, are doubled as compared to the case of Mi ~ m, 
i.e. we get 



Kir) - -Clo,-1- 



(37) 



for M2 — > 00. 

The only essential result of changing AI2 from m to infinity is that the critical elastic constant increases by a factor 
of two. 



IV. SUMMARY AND CONCLUSIONS 



For a translationally invariant chain with one anharmonic bond and otherwise harmonic nearest neighbor interac- 
tions, we have shown exactly how the anharmonic degree of freedom can be separated from the harmonic ones in their 
normal mode representation. As a result, we have obtained Lagrangian (21), which is of the form of Lagrangian (4). 
Note, that this result can also be obtained for a three-dimensional system within the Born-Oppenheimer approxima- 
tion starting with an arbitrary translationally invariant potential V{xi, . . . ,xn) for a iV-particle system jl2l |. Since 
the Caldcira-Leggett type nonlocal action [U, [3, Q is based on the form (4) (or (21)) of the Lagrangian it is not the 
translation invariance and therefore not the conservation of momentum, which can lead to a different type of nonlocal 
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action. The discrepancy between our results and those of Ref. [3| niay have the following origin. Since the harmonic 
part of Lagrangian ([7]) is diagonal in Xa, the harmonic variables Xa are already normal mode coordinates. In that 
case a translation of the full system only changes the Goldstonc mode amplitude, let us say xq, but leaves all the 
other normal mode coordinates unchanged, ^ for any translation. If xi and X2 in Eq. ([7]) are real 

space coordinates then the coupling term {xa — X2)^ is not translationally invariant for a 0. 

Our model has allowed us to calculate explicitly, e.g. for m„ = m and AI = N/2 the coupling constants Ca, the 
eigenfrcqucncics cUa and the spectral density J{uj). The w-dcpcndencc of J is determined by the density of states g{ijj) 
and the coupling constants Cq. Although for a; is independent of the number M, the frequency dependence 
of Cq, exhibits a sensitivity to the location M, which makes the affect of the harmonic bath on quantum tunneling M 
sensitive. As a consequence, the damping is Ohmic if the bond defect is within the bulk of the chain and superohmic 
if it is close to the boundaries. For the former case there is a transition from a delocalized state (due to tunneling) to 
a localized one if the elastic constant exceeds a critical value Ccrit, whereas tunneling is never suppressed in the latter 
case, provided the 'observation time' T is large enough compared to tm which is roughly M times the instanton kink 
width. For T < tm (since the thermodynamic limit N ^ oo had already been performed, M must be finite, but can 
be arbitrary large) the dissipation is effectively Ohmic leading to a similar behavior when the bond defect is within 
the bulk. 

If M ~ 0{N) (e.g. M = N/2) and if one of the masses of the bond defect tends to infinity, no significant changes 
occur except for doubling of the critical constant Cent- This is obvious since, e.g. M2 — ^ 00, makes the part of the 
chain to the right of the defect inactive, i.e. the phonons to the right do not act as a reservoir for the bond defect. 
Accordingly, only half of the harmonic chain is generating dissipation, which results in doubling of Ccrit ■ 

Although we are not aware of a concrete experimental system, these results could be relevant for linear macro- 
molecules, which may be described by the model Hamiltonian Eq. (10). If many of such molecules with a single defect 
are produced, the position of which can be controlled experimentally, one might observe, e.g. the location-dependent 
tunneling by spectroscopic methods. 
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APPENDIX A: USE OF COM AND RELATIVE COORDINATES OF THE TOTAL CHAIN: 

DIAGONALIZATION OF THE MATRIX T 



The separation of the harmonic and anharmonic DOF by using the center of mass and relative coordinates of the 
total chain has been described in Section [Til The transformation to normal coordinates requires the diagonalization 
of the matrix T. Here the most important steps of the diagonalization procedure are outlined. 

Making use of Eqs. (12)-(15) one obtains the matrix elements of the symmetric matrix T in the form 



and 



( mi + nii+i 
miini+i 

mM-i + {mM + mM+i) 
mM-i{mM + mM+i) 

{■niM + niM+i) + mM+2 



l,...,A/~2,M + 2, ...,iV-l, 



i = M - 1, 




M + h 



l,...,M-2,M-M,...,iV-2, 



i = M -1 



i+2 



1 



mM + mM+i 



l,...,M-3,A'/-|-l,...,iV-3, 
M -1 . 



(Ala) 



(Alb) 



(Ale) 



12 



The diagonalization of T can not be done analytically for arbitrary masses m^. Therefore we take the simplest case 
of equal masses, = m. Then it is straightforward to prove that the eigenvalue equation 



N-l 



is solved by 

{sin(g„i), I < i < M - 1, 

(A2) 
baSm{qaiN - i)), M + l<i<N-l, 

A„ = -(l-cos(g„)) (A3) 
m 

where Afa is the normalization constant and ba is a coefficient depending on the location AI of the bond defect. The 
wave numbers qa are solutions of the transcendental equation 

cot(gjV) = cot(gA/) + . . (A4) 

2 sm [qAl ) 

Since the l.h.s. of Eq. (jA4|) diverges at g = • a, a = 0, 1, ... it is easy to see that its solutions are of the form 

TT 

qa = —a + ea, a = 0, 1, - 3 (A5) 

with < ea < jj- There are (A^ — 2) solutions corresponding the (A^ — 2) harmonic DOF. The remaining two DOF 
are the COM and the bond defect coordinate Xc and qM, respectively. In the thermodynamic limit A^ — *■ oo, the 
discrete set of qa wave vectors becomes a continuous variable q within [0, tt] with a constant density which together 
with Eq. (jA3P implies a constant low energy density of states. 

The normalization constant Afa and the coefficient ba are functions of qa, AI and A^. Their explicit expressions are 
not given here. 

APPENDIX B: USE OF COM AND RELATIVE COORDINATES OF THE BOND DEFECT 

In this appendix we will describe the separation of harmonic and anharmonic DOF using an approach alternative to 
that used in Section |TT1 It has the advantage that it can be straightforwardly applied to higher dimensional systems. 
The starting point is the introduction of COM and relative coordinate and gM, respectively, of the bond defect 
(see Eqs. (|lla|) and (jllb|l ). The corresponding canonical momenta 

Pd = Pai + PM+i, (Bla) 

■KM = ; Pm+1 ; Pm- (Bib) 

niM + mM+i ruM + rnu+i 

Substituting Xd, qu from Eq. (11) and Pd, ttm from (Bl) into Eq. (10) yields 

H = Hd + Hharm + Hint (B2a) 

where 

Hd = + yoilM) + -777 72lM B2b) 

1 , 1 , c 



2m„ 2{mM + mM+i) 2 ^ 

(n^J\/,J\/ + l) (n^A/,J\/±l) 



13 



C C 



(B2c) 



Hit 



-C 



rriM+i 



rriM + niM+i 



{Xd — xm-1 — flAZ-i) + 



TOM 



TO A/ + TOM+1 



-{xm+2 — Xd~ flM+l) 



qM 



(B2d) 



Here /j,m = rnn'mM+i/ {mM + wa/+i) is the reduced mass of bond defect. Note that Hd from Eq. (jB2bp is identical 
to Hd from Eq. (jl6bp after replacing (ttm, (Zm) by (pm, 9a/)- The transformation of Hharm (Eq. (|B2cP ) to the normal 
coordinates may be more conveniently carried out using the notations 



Xn , n = l,...,M-l 

Xd , n = M 

Xn+i , n = AI +1,...,N -1 



(B3a) 




n = 1,...,M -1 

n = M 

n = M +l,...,N - 1 



(B3b) 



and 



i„ , n = l,...,Af-l 

rrLM + raM+i , n = M 
TUn+l , n = il/ + 1, iV - 1 

Next WG expand the potential part \^harm{'^i7 ■■■;-^7v- 

of Hharm around its equilibrium configuration, 



(B3c) 



(B4) 



up to the second order terms in it^. Note that this is not an approximation since Vharm is a harmonic potential. This 
leads to 



JV-2 



Hharm = ^ 2" ^ ^""+^ ~ ""^^ 



(B5) 



Introducing the mass-weighted coordinates 



(B6a) 



(B6b) 



Eq. dH]) yields 



7V-1 



1 ''"^ 1 ^"^ 

Hharm ~ ^ Pii + ^ ^nrn^n^rn 



2 2 

n— 1 n.rn— 1 



(B7a) 



where the only nonzero matrix elements of the symmetric matrix V are 



' r 



f 

nn 



C I 1 , n = l,iV- 1 
<|2 , n = 2,...,iV-2 



(B7b) 



and 



* nn+1 



C 



(B7c) 
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Let el""* and A^, a = 0, 1, TV — 2 be respectively the eigenvectors and eigenvalues of V. The canonical transfor- 
mation 



N-l N-1 

= ^ ^ '^n^n \ Pa — ^ ^ Pn^n ^ 
n=l 



n=l 



leads to the normal mode representation 



N-2 



harm 



a=0 



(B8) 



(B9) 



Since Hharm in Eq. IjBSp is still translation invariant there is a zero frequency mode which we choose for a = 0. With 
Aq = we get 



N-2 



Hharm = ^Po + ^ 



+ AqX^ 



(BIO) 



The first term in the r.h.s. of Eq. (jBlOp is the kinetic energy of the COM of total chain. The second term corresponds 
to Hharm from Eq. mil). Using Eqs. (B3), (HU, (B6) and dHI brings the interaction part (Eq. (jB2dp ) into the 
form 



Hi, 



(Bll) 



with 



C- 



rriM + ruM+i 



fTlM+l 



-.(a) 



M 



~(a) 




M+1 /m~ M 



M 



(B12) 



Again, the analytical diagonalization of V cannot be performed for arbitrary masses. Accordingly, we choose m„ = m 
as in Appendix El Then Eqs. (jB7bp and (jB7cj) result in 



II , n^M, 



n = M, 



(B13) 



m 2 



n= l,il/,iV-- 1, 
n^^ l,A/,iV- 1, 



(B14a) 



m l/v2 



n^M - 1, Af, 
n = M - 1,M 



(B14b) 



All the other matrix elements vanish. Again it is straightforward to prove that the eigenvalue equation 
El=l K,„e^^^ = Aaei"^ for n^M- 1, Af, M + 1 is solved by 



g(")=A/' /cos (9a (n- 5)) , n = l,...,Af-2, 

" UaCOs(qa (iV-n- i)) , 71 = A/ + 2, ...,7V- 1, 



(B15) 



Aq = —(1 - cos(gc<), (B16) 

with J\fa being the normalization constant and ba a A/-dependent coefficient. The remaining equations for e^"'' with 
n = Af — 1, Af and Af + 1 yield a nontrivial solution if a corresponding determinant vanishes. This condition leads to 
the transcendental equation 



(B17) 
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for the wave numbers g^. Although Eq. (jB17p looks quite different from the transcendental equation (jA4[) it can be 
shown by use of identities for trigonometric functions that (|B17p and (|A4p are equivalent, i.e. the set of solutions 
{qa} of Eq. (jB17[) and {qa} of Eq. (jA4p are identical. We have already stressed that Hd from Eq. (|B2bp and that 
from Eq. (|16bp are identical, as well. Straightforward but tedious calculations show that the complete Lagrangian 
corresponding to the classical Hamiltonian Eq. (B2) is identical to the Lagrangian (21). Particularly, it can be proven 
that Cq from Eq. (|B12p is identical to Ca from Eq. ([20|) . 



APPENDIX C: SEPARATING THE HARMONIC PART INTO LEFT AND RIGHT PARTS 



In this appendix we will show that separation of the harmonic and anharmonic DOF can be done by taking the 
left and right harmonic parts separately. Similarly to the approach used in Section [ill our first step is to separate the 
COM of the total chain from the relative coordinates. This leads to the Hamiltonian from Eq. Neglecting the 

kinetic energy of COM Eq. ((T4|) can be rewritten as 

H ^ Hd+ Hf^^„^^ + H^arin + Hint (Cla) 

where 

Hd = -^'Kll + V^{qM), (Clb) 



M-1 „ M-1 

tL ^ \ " rj^L , ^ \ " „2 



HLrn. = 2 Tt^^^^j +-^qf^ (Clc) 



iV-1 ^ N-1 

ij=A/+l " i=M+l 



^ N-1 „ Af-1 



Hint = — ( ttm-i H ttaz+i ) ttm (Clc) 

\mM rriM+i J 

and the nonzero matrix elements are 

J-ii - : ^ij+i ^i+i.ti o-L,K (Cltj 

with i = 1, A/ - 1 for cr = L and j = 7\/ + 1, TV - 1 for cr = i?. Let ef ^''^ (ef ^^^) and A^(A^) be the eigenvectors 
and eigenvalues of T^(T^). 
Then we use the notations 



X 



M-1 N-l 



i=l i=M+l 



M-1 N-l 
i=l i=M+l 



in order to get 



N-l 



fj,=M+l 
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for the harmonic part and 



Hi, 



' A/-1 
i/=l /i=M+l 



N-1 



(C4a) 



for the interaction with 



rriM 



mM+1 



-A/+1 



(C4b) 



for the couphng constants. 

This type of approach describes the chain as a bond defect coupled to two baths of harmonic oscillators, the left 
and right part of the chain. For the path integral formalism we need the Lagrangian. From (|Cla[) . (|Clb|) . (C3) and 
(|C4aP we can determine the velocities qM,ii and as function of the momenta. Solving for the momenta as a 
function of the velocities is straightforward but tedious. We report the final result 



■KM = K 



M-1 r N-M-1 R 



qM 



(C5a) 



= —x^ 



N-M-l 

=1 ^M' 



E 



(C5b) 



1 



A/-1 r N-M-l R 



(C5c) 



with 



1 - MA/ ^ 



A^ 



E 



(C6) 



Making use of (C5) for the calculation of the Legendre transform of H from Eq. (CI) leads to the Euclidean Lagrangian 

L = Lci + Lfiarm + -^int (C7a) 

where 

Ld^^ql + VoiqM) (C7b) 

in which 



A/-1 



1 

A^ 



7V-A/-1 



A/-1 

E 



A^- 



Af-A/-1 p • 



(C7c) 



JV-A/-1 ^fl 



AZ-l ^ 



(C7d) 



This form of L differs completely from that of Eq. (21). Particularly, Li^t from Eq. (jC7dp is a coupling of the 
velocities and not of the bond defect coordinate qM with the normal mode coordinates Xa as for Lint from Eqs. (|21ap . 
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(|21c[) . In addition, the harmonic part Eq. (|C7c|) is not 'diagonal', i.e. due to the third term on the r.h.s. of Eq. (|C7c|) 
there is an intra- and an inter-couphng between the phonons (normal modes) of the left and right harmonic part of 
the chain. 

In order to eliminate the harmonic degrees of freedom in the path integral representation of the propagator one has 
to 'diagonalize' Lharm from Eq. (|C7c[) . This can be done by a point transformation x{^{{xa}, Qm) and xj}{{xa}, Qm)- 
This transformation follows directly from Eqs. ()15bp . (18) and (C2): 

N~2 /M~l \ 

xi{{xa},qAi) =J2\Y1 Sj^^ef^"^ Xa ^HlM+l e^|'l\qA/ (C8a) 

Q— 1 \ %—l / 



N-2 / JV-1 \ 

J, m) = E E ^'^f''' - ^ 71 ^ZWf (csb) 

a=l \i=M+l / 

Taking the time derivative of Eq. (C8) yields the transformation of the velocities. Substituting this and the trans- 
formation (C8) into Eq. (C7) 'diagonalizes' Lharm and replaces the velocity coupling by a coupling of qM and {xq}. 
After a lengthy calculation one arrives at the Lagrangian from Eq. (21), which, of course, is not a surprise. 
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